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Abstract 

We present an extension of the model proposed recently to account for dramatic 
changes below T c (anomalies) of some c-axis polarized infrared-active phonons in 
bilayer cuprate superconductors, that applies to trilayer high-T c compounds. We dis- 
cuss several types of phonon anomalies that can occur in these systems and demon- 
strate that our model is capable of explaining the spectral changes occurring upon 
entering the superconducting state in the trilayer compound Ti2Ba2Ca2Cu30io. 
The low-temperature spectra of this compound obtained by Zetterer and coworkers 
display an additional broad absorption band, similar to the one observed in un- 
derdoped YBa2Cu307_<5 and t^S^CaC^Os- In addition, three phonon modes are 
strongly anomalous. We attribute the absorption band to the transverse Josephson 
plasma resonance, similar to that of the bilayer compounds. The phonon anomalies 
are shown to result from a modification of the local fields induced by the formation 
of the resonance. The spectral changes in Tl2Ba2Ca2Cu30io are compared with 
those occurring in Bi2Sr2Ca2Cu30io, reported recently by Boris and coworkers. 
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1 Introduction 



Some of the phonon modes in the high-T c cuprate superconductors are strongly 
renormalized when going from the normal to the superconducting state [1,2,3,4,5,6,7,8,9,10]. 
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It is important to establish what particular type of electron-phonon interac- 
tion is responsible for these effects (so called phonon anomalies), and whether 
the anomalous phonons are "active players" in the mechanism of the high-T c 
superconductivity [11,12,13] or whether they simply respond to changes of the 
electronic ground state that have an independent cause. Here we focus on the 
anomalies of c-axis polarized infrared-active phonons. The most pronounced of 
them have been observed in the so called multilayer compounds having two or 
more copper-oxygen planes per unit cell. Recently Munzar et al. [14] explained 
phonon anomalies in underdoped YBa 2 Cu 3 07_5 (Y-123). The explanation is 
based on the multilayer model of van der Marel and Tsvetkov [15,16,17], and 
the assumption of a week (Josephson) coupling between the copper-oxygen 
planes: it has been assumed that the interlayer conductivities are fairly inco- 
herent in the normal state and dominated by the coherent superfTuid contribu- 
tion in the superconducting state. A bilayer superconductor is thus considered 
as a superlattice of inter-bilayer and intra-bilayer Josephson junctions — this 
picture is called the Josephson superlattice model. The presence of two dif- 
ferent junctions in a unit cell gives rise to a new transverse mode, the so 
called transverse Josephson plasma resonance (t-JPR) with a frequency be- 
tween the one of the inter-bilayer (longitudinal) and the intra-bilayer (longi- 
tudinal) plasmon. The same picture of the c-axis charge dynamics of a bilayer 
system has also been arrived at by N. Shah and A.J. Millis [18] starting from 
a microscopic theory involving in-plane Green functions. An important new 
ingredient of the model of Ref. [14] is that the local field effects are taken 
into account: the changes of the phonon resonances in the infrared spectra 
result from changes of the (dynamical) local electric fields associated with the 
formation of the t-JPR. The same model, with a slightly different parameter- 
ization of the electronic contribution, has been more recently used to explain 
phonon anomalies in another bilayer compound, Bi2Sr2CaCu20 2 (Bi-2212) [5]. 
Very recently the t-JPR and the related phonon anomalies have also been ob- 
served in the trilayer compound Bi 2 Sr 2 Ca2Cu30io [6] and briefly discussed in 
terms of the Josephson superlattice model and the local field effects. Here we 
present a detailed account of several types of phonon anomalies that can occur 
in trilayer cuprate superconductors (Sec. 2) and we compare our predictions 
with experimental data for the trilayer compound Tl 2 Ba 2 Ca2Cu30io(Tl-2223) 
obtained by Zetterer et al. [7] (Sec. 3). Surprisingly, these data are found to 
be fairly different from those of Bi 2 Sr 2 Ca2Cu30io. In Sec. 4 we summarize our 
conclusions. 
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2 Theory 



2.1 The basic concept and the response of a stack of superconducting CuOi 
planes 



A schematic representation of the three Cu02 planes of a trilayer compound 
and the basic structural element of Tl-2223 are shown in Fig. 1; E h \ and E int 
denote the average electric fields in the intra-trilayer and the inter-trilayer re- 
gion, respectively, and k (—k) denotes the surface charge density of the lower 
(upper) outer CuC>2 plane. The distance between the copper-oxygen planes in 
the trilayer block is denoted by d^\ and the distance between neighbouring tri- 
layers by d int . The magnitude of the inter-trilayer current density j int is, due to 
the low conductivity of this region, much smaller than that of the intra-trilayer 
current density, jbi- Consequently, ji nt will be neglected (ji nt = 0). Further- 
more, the current density jbi can be expressed as = —iujeoXb\{^)Eu{^), 
where Xw(^) is the local susceptibility of the Josephson junction between the 
closely spaced Cu0 2 layers. Because of symmetry, the local susceptibility is the 
same for both regions inside the trilayer block. As a consequence, the middle 
cooper-oxygen plane does not become charged in an infrared wave; the fields 
and the current densities in the two regions are the same. The trilayer block 
behaves as one Josephson junction with the susceptibility Xbi(^) connecting 
the outer planes of the block. 

The c-axis dielectric function in the optical limit can be written as (cf. [19] p. 
127) 

e(w ) = i + -LiM, (1 ) 

6qUJ E(0l>) 



where and E{uj) are Fourier components of the total induced current 

density along the c-axis and the total electric field, respectively. Both vari- 
ables are macroscopic, i.e., unit cell averages of the corresponding microscopic 
quantities. This equation can be written as 

«M = «- + , (2) 



where is the interband dielectric constant, the induced current density is 
expressed as the sum of microscopic current densities ]\{uj) and the volume 
average is explicitly denoted by (). For a trilayer compound, Eq. (2) can be 
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further rewritten as 

€{U}) 600 + 24i + dint e Q u E(u) + e c E{u) ' 1 J 



where the term 2<ibi/(2<ibi + dint) is the volume fraction of the intra-trilayer 
region and j% h (ou) = —iumkQk r k(u) represents the (volume averaged) current 
density due to vibrations of the ions of type k. Further n&, Qk and rj,(uj) are the 
corresponding concentration, effective dynamical charge, and displacement, 
respectively. 

First we discuss the electronic part of the model, i.e., we set j! h = 0. The 
electric fields E b \ and E int shown in Fig. 1 a) are determined by the following 
equations 

E hl (u) = E'(uj) + , (4) 

e e oo 

E int (iu)=E'(Lu) , (5) 
EH = -^—E^H + -Jf—E^H , (6) 

ZCtbl + «int ^«bl + "int 



where E' is the component of the average internal field E that is due to external 
charges and depolarization processes at high frequencies. In the absence of the 
latter (i.e. for eoo = 1) E' would correspond to the field Eq of Kittel's text 
book [20]. The continuity equation reads 

j u (u>) = \ujh(uj) . (7) 



Equations (4)-(7) are solved by 

E iQt {uj) (2d u + di n t) ebi(^) E u (u) (2d u + d int ) e c 



(8) 

/ m\UJ) ^u-bltoo T "inttblV^; 

. (24l + dint) M 



E(u) 2dbi£t» + d int e b i(cj) ' E(u) 2d b ie 00 + d- mt e hl (uj) 

[2d u + di n t) e bl 
2cibieoo + dintebi(^) 



where e M (w) = + XuM- 

In the following, the susceptibility xw is taken in the form Xw( w ) — ~ UJ b\/ uj2 + 
iSbi/u. The first and the second term represent the response of the super- 
fluid and the background, respectively 1 . The background is described in the 

1 It has been recently suggested [21] that the first term should be given a finite 
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simplest possible way (an infinitely broad Drude term) in order to keep the 
number of parameters small. Figure 2 shows the spectra of the real part <7i(u;) 
of the optical conductivity cr(ui) = — icje e(cj) for several values of the plasma 
frequency u>u- Since this article aims at explaining the experimental data for 
Tl-2223, the following values characteristic for this compound have been used: 
du = 3.2 A, di n t = H-5 A, A = a 2 , where a = 3.9 A is the in-plane lattice 
parameter [22]; further eoo = 5 (the value typical for the c-axis response of 
the high-T c cuprates), and = 1400 cm' 1 (a value, that gives a reasonable 
width to the t-JPR). The maximum in the spectra of Fig. 2 corresponds to 
the t-JPR. It can be seen that a^i determines not only the frequency of the 
resonance but also its spectral weight. 

In order to include phonons, we have to take into account the local fields acting 
on the ions participating in the vibrations. It can be seen from Fig. 1 a) that 
the charge densities k and — k give rise to three different electric fields: (a) the 
field between the trilayer blocks (E\°^), (b) the one inside the trilayer (E^ c ), 
and (c) the field at the outer copper-oxygen planes (E^£L). They are given as 



As a consequence, we arrive at three different elementary types of phonons 
consisting in vibrations of ions located in the three regions. 



2.2 Vibrations of ions located in the inter- or intra-trilayer region 

In order to include a phonon mode consisting in a vibration of a given type 
of ion, we have to express its contribution to the current density, j ph , and 
describe its influence on the fields Eu and E int - For concreteness we focus on 
the inter-trilayer case; the intra-trilayer one can be treated in an analogical 



spectral width, since it could be unrelated to the SC condensate. The suggestion 
was motivated by the observation that in strongly underdoped Y-123 the additional 
resonance, which is determined by this term, appears already at temperatures much 
higher than T c . In other compounds, however, it starts to appear below T c or very 
close to T c . We are convinced that at least in the superconducting state the first term 
corresponds to the condensate. The origin of the above- T c phenomena in strongly 
underdoped Y-123 is not yet clear but we believe that they are due to fluctuation 
effects. 




(11) 



(10) 
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way. The classical equation of motion for a particle of mass m and charge Q 
yields the following formula for the Fourier component of the displacement 



r{uj) = 



nQ 



(13) 



where x is the polarizability of the phonon, 

X{u) = — Wlth S Q = 2 ■ 14 ) 



Further n is the density of the vibrating ions, n = N/[A(2d h \ + d int )}, where N 
is the number of the ions per volume A{2d^\ + dint); u q an d 1q are the char- 
acteristic frequency and the broadening parameter, respectively. The (volume 
averaged) current density corresponding to this vibration is 

f\u) = -kjeo X (u) E£(u) . (15) 



The influence of the phonon on the electric fields can be taken into account 
using the following simple model: the planes of the ions are approximated by 
uniformly charged planes perpendicular to the c-axis with the same surface 
charge density, rj = NQ/A (see Fig. 3). Since the electric field generated by 
a uniformly charged plane does not depend on the distance from the plane, 
this assumption simplifies the model considerably. Using Fig. 3, we obtain the 
following equation for the average field in the inter-trilayer region: 

4rt£int(f) = K + r(t)](-E v ) + [di^ -dr,- r(t)]£„ + d mt E'(t) , (16) 



where d v is the equilibrium distance of the charged plane from the lower CuC>2 
plane, and E v = r]/2eoeoo- Screening by high frequency processes is taken into 
account by in the denominator. Fourier transformation of Eq. (16) yields 



E- mt (uj) = E'(uj) 



r{uj) rj 



dint ^0 Coo 

7 / x cloc/ ,\ . 2rf bl + dint 



^H-^xHW, 7= P ' mt - (it) 

For the particular type of phonon Eq. (17) substitutes for Eq. (5), whereas 
Eq. (4) remains unchanged. The second term on the right-hand side of Eq. (17) 
represents the depolarization field of the phonon. Note that it contributes 
only to the field E int and not to the field E u . On the other hand, the volume 
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averaged depolarization field 

"XM^M (18) 



rf int 7 . . , w ,lo< , 



2du + d- m t 



is equal to — P/e eoo ; where P is the macroscopic polarization [P = e x(c<j)£ , i 1 ° t c ]. 
This is the textbook formula for the macroscopic depolarization field {E\ of 
Kittel's book) of a thin plate. The local fields E\°%, E^ c , E 1 ™ do not contain 
any contribution of the depolarization fields. This resembles the exact cancel- 
lation of the depolarization field and the Lorentz field (E 2 of Kittel's textbook) 
that occurs in cubic crystals. Note that the phonon polarizabilities x entering 
the model equations represent response functions with respect to the local 
fields instead of the averaged field. In the absence of interlayer currents, the 
input frequencies would correspond to the LO-frequencies while the frequen- 
cies renormalized according to the model equations would correspond to the 
TO-ones; considering only one phonon we obtain uto = — S^/e^. Note 

that if we include the current j ph in Eq. (3), we also have to take into account 
the influence of the corresponding polarization on the electric field, i.e., the 
depolarization field. Otherwise the approach would not be consistent and may 
lead to unphysical results. 

Figures 4 a) and 4 b) show the real part of the optical conductivity for two val- 
ues of the phonon frequency, ujq = 320 cm -1 and uq = 610 cm -1 , respectively. 
The values of the remaining parameters correspond to 02 or 03 oxygens in 
Tl-2223: N = 2, m = 16amu, and we take Q = — 2|e| and 7q = 20 cm -1 . 
The full lines represent the spectra of the normal state (NS) with cou = and 
display only the phonon resonance. The dashed lines represent the spectra of 
the superconducting state (SCS) with uu = 1400 cm -1 and display both the 
t-JPR (the additional peak) and the phonon. We emphasize that the set of 
the values of the parameters used to compute the SCS spectrum differs from 
that of the NS only in the value of the parameter uu- It can be seen that if 
the frequency of the phonon is lower than that of the t-JPR [Fig 4 a)], the 
spectral weight (SW) of the phonon is enhanced in the SCS. If it is higher 
[Fig 4 b)], the SW is reduced. Note that the qualitative characteristics of the 
anomaly depend only on the order of the two frequencies. 

Next we consider a phonon consisting in a vibration of intra-trilayer ions of a 
given type. The derivation of the formula for the depolarization field is entirely 
analogical to the case of the inter-trilayer ions. We obtain 



E u (uj) = E'(lo) + 
= E\uj) + 



e e oo 
£0 e oo 



2(4,1 eo e oo 
5 



2gLi + di 



int 



2d 



bi 



(19) 
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This equation replaces Eq. (4) whereas Eq. (5) remains unchanged. The main 
difference from the inter-trilayer case is that the ions feel the local field E]^ c 
instead of E}°^. Figures 4 c) and 4d) show, how the phonon resonance is influ- 
enced by the superconducting transition. The values of the phonon frequencies 
are: ujq = 320 cm -1 and 610 cm" 1 , respectively. The values of the other param- 
eters correspond to Ca in Tl-2223: N = 2, m = 40amu, and we take Q = 2|e| 
and 7q = 3 cm -1 . It can be seen that if the phonon frequency is lower than 
that of the t-JPR [Fig 4 c)], its SW is reduced and the phonon softens consid- 
erably. If it is higher [Fig 4d)], the SW is enhanced and the phonon frequency 
increases. Note that the phonons are renormalized already in the NS because 
of the background term in the formula for xu- The renormalization manifests 
itself, among others, in a non-Lorentzian shape of the resonances. 

We have shown in Sec. 2.1 that the trilayer block behaves simply as a bilayer 
one with the interplanar distance of 2<i bl . The approach presented here can 
thus be applied to both trilayer and bilayer compounds: in the latter case, 
2dbi in the coefficients 7 and 5 has to be substituted by d\,\ and we have 
7 = (4i + d int )/d int and 5 = (d hl + d int )/d u . 



2.3 Vibrations of ions located in the outer CuC>2 planes of the trilayer block 



For concreteness we further discuss the infrared- active (in-phase) vibrations 
of the oxygens 01 of the outer Cu02 planes. Vibrations of the copper ions, 
however, can be treated in the same way. Boundaries between the intra- and 
inter-trilayer regions can be identified either with the static planes of the Cu 
ions or with the planes of the vibrating oxygens. In order to simplify our 
considerations we choose here the first approach. We have checked that for 
reasonable values of effective charges the use of the second one leads to very 
similar results. When expressing the depolarization fields in any of these ap- 
proaches we encounter a problem: the interface between the inter- and the 
intra-trilayer regions is repeatedly crossed by a charged plane, which makes 
the time dependence of the fields and E int unharmonic. This difficulty 
results from our approximation of the ionic planes by two-dimensional homo- 
geneously charged planes. In order to get rid of the problem and keep the 
useful approximation at the same time, we represent the planar oxygens 01 
by two homogeneously charged planes as shown in Fig. 5. 

The surface charge density of each vibrating plane is rj/2 with rj = N01Q01/A, 
Noi = 2. Following the derivation of the preceding section, we obtain the 
subsequent formula for the average field between the upper two Cu02 planes 
of the trilayer: 
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BblM = £V) + ^-^^- = 

eo^oo «bi eo e oo 

= s'H + xoiM ^ocp(^) , « = — T-i — • (20) 

Co^oo Coo 4(3bl 

We have used the relation between the displacement r and the local field 
E l °c P , r = eoXoi-^ocpA^oiQoi) [ c ^- Eq. (13)], where n i is the concentration 
of the oxygens 01 of an upper Cu0 2 plane, n i = N i/[A(2d h i + d int )}. The 
average field between the two lower Cu0 2 planes is also given by Eq. (20). 
Consequently, the central Cu02 plane remains without time dependent charge 
density even for this type of vibration. The equation for the field E int , 



dint £0 e c 

= ^xoiM^M, (3= 2 -^±^ (21) 



can be obtained in an analogical way. Equations (20) and (21) substitute 
for (4) and (5). The current density due to the vibrating oxygens is given as 
the sum of the contributions of the two outer oxygen planes [cf. Eq. (15)] 

f\uj) = -2\ue oX o^)E^{u) . (22) 



The formulas for a and (3 of a bilayer system can be obtained along the lines of 
the considerations at the end of the previous subsection: a = (rfbi + di nt )/2d h \ 
and = (4i + d mt )/2d- mt . 

Note that the second approach (with the boundaries associated with the 
charged planes of the planar oxygens) yields slightly different expressions for 
a and (3 in Eqs. (20) and (21) and also the formula for the current density 
due to the phonon differs from Eq. (22), see Appendix. In the previous works 
concerning bilayer systems Y-123 [14] and Bi-2212 [5,23], a and f3 have been 
calculated using the second approach but j ph using the first one, i.e., Eq. (22). 
Fortunately, the errors caused by this formal inconsistency are negligible be- 
cause the values of a and j3 resulting from both approaches are very close. 
Note that there are some errors in the third line of Table 1 of Ref. [14], the 
correct values will be published in an erratum [24]. 

Figure 6 a) and 6 b) shows the interplay between the t-JPR and the phonon 
for two values of the phonon frequency: u; i = 400 cm" 1 and 660 cm -1 . The 
values of the remaining parameters are iVoi = 2, m G i = 16amu, and we 
take Qo\ = — 2|e| and 701 = 20 cm" 1 . If the phonon frequency is lower than 
that of the t-JPR [Fig 6 a)], the formation of the t-JPR is associated with a 
reduction of the phonon SW and its softening by about 5 cm -1 . An even more 
pronounced phonon anomaly of this type occurs for the bond bending mode 
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in underdoped Y-123 [14]. The fact, that the anomaly expected to occur in 
Tl-2223 is smaller than in Y-123, is connected to the difference in the width of 
the multilayer block. If the phonon frequency is higher than that of the t-JPR 
[Fig 6 b)], the SW of the phonon is enhanced in the SCS. 



2.4 Composite vibrations 



Up to this point we have described elementary vibrations where only one 
type of ion participates. This can be a good approximation for a mode where 
the major part of polarization is due to vibrations of one type of ion. For 
modes where several types of ions contribute significantly an extended version 
of the model is required. In this section we study composite modes of the 
oxygens of copper-oxygen planes. According to the results of the shell model 
calculations [26] , in trilayer compounds such modes can be expected to occur. 

In order to treat the composite modes as simply as possible we introduce the 
mechanical coupling scheme shown in Fig. 7, which provides an out-of-phase 
and an in-phase resonance. The equations of motion read 



d 2 r(t) dr(t) 
m o^r- = -hr(t) + k 2 [z(t) - r(t)} - m oloi ^- + F r (t) , (23) 

m °^d^ = 2fe[r(t) ~ ^ ~ mo7 ° 4 ^ + Fzit) ■ (24) 

The forces F r and F z are proportional to the local fields: F r = E^pQoi and 
F z = E^fQoi. For the Fourier components of the ionic displacements we 
obtain 

F r (u) q z {uj) + F z (u) k 2 F z {uj) q r {uj) + 2F r (uj) k 2 /n _, 

r y jj > = r~\ — r~\ — oTi — > z y u ) = r~\ — r~\ — 072 — • I 25 ) 



Here q r = —u) 2 mo + k\ + k 2 — \uJ^ioi m o > qz = —0J 2 mo + 2k 2 — icJ7o4^o- 
The depolarization field for a composite vibration is the sum of the contribu- 
tions of the elementary vibrations. Following the consideration of the previous 
subsections we obtain the equations for the mean fields: 



E h M = E'{uj) + rM*7Qi + *Mqo4 (26) 

^o e oo 2abieo e oo 

BlntM= ^-lM*L, (27) 

«int ^0 e oo 

where 7701 = NqiQqi/A and 7704 = Nq^Qq^/A. Equations (26) and (27) again 
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replace Eqs. (4) and (5). The current density due to the vibrations of the 
oxygens is 

jfc h (^) = -2iun 01 Q ir(uj) - wn AQoAz{uj) ■ (28) 

k 

The first and the second term on the right side of the Eq. (28) is due to vi- 
brations of oxygens in the outer and inner copper-oxygen planes, respectively 

Figure 8 shows the real part of the optical conductivity for the model described 
above. The values of the parameters are: N Q i = N 0i = 2, Q 01 = Q Qi = 
— 2 1 e | , mo = 16, and we have chosen 701 = 704 = 20 cm -1 , k<ijk\ = 0.3 and 
Ul = = 320cm- 1 (a), 600cm- 1 (b), 800cm- 1 (c), and 1100cm- 1 (d). 
The model yields an out-of phase (OPR) resonance and an in-phase (IPR) 
one. If the OPR frequency is lower than that of the t-JPR [Fig. 8 a)], its 
SW increases considerably when going from the normal to the SCS. On the 
contrary, if its frequency is higher [Fig. 8 b)], the SW of the phonon decreases. 
In both cases the phonon frequency remains approximately constant under 
the superconducting transition. For the IPR the situation is different: if its 
frequency is lower than that of the t-JPR [Fig. 8 b) and c)], its SW is strongly 
reduced and its frequency decreases. If the IPR frequency is higher [Fig. 8d)], 
the changes are reverse: the SW increases and the phonon hardens. In Fig. 8 a) 
the interplay of the depolarization fields and the electronic background leads 
to a strong overdamping of the IPR. In Figs. 8 c) and d), the OPR is located 
above the displayed frequency range. It can be seen that the IPR frequency 
is renormalized much more than that of the OPR. This can be understood: 
for the IPR the displacements of the oxygens r and z have the same sign, 
therefore the depolarization field (the third term on the right-hand side of 
Eq. (26), determining the frequency shift, is large. On the contrary, for the 
OPR the signs of the displacements are opposite and the depolarization field 
is therefore small. Note that the qualitative aspects of the phonon anomalies 
of the composite modes are insensitive to changes of the value of the ratio 
h/k 2 . 

The changes of the spectral weights reported in this section can be elucidated 
in terms of the changes of the local fields. Figure 9 displays frequency depen- 
dencies of the ratios E^ c /E, E\°^/E and E l ° cp /E in the SCS obtained using 
Eqs. (4)-(12). The phonons have not been taken into account and Su = 0. 
The singularity at u p ~ 500 cm -1 corresponds to the t-JPR. Since the values 
of the ratios in the NS are equal to one, their values in the SCS express the 
relative changes of the local fields when going from the NS to the SCS. Let 
us consider, e.g., the local field E\°^ that acts on ions situated in the inter- 
trilayer region. If the frequency uoq of a phonon involving vibrations of these 
ions is lower than u p , the local field E\°^(uq) increases and the SW of the 
phonon increases [cf. Fig 4 a)]. If the frequency is higher, EI°^(uq) decreases 
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and so does the spectral weight [cf. Fig 4 b)]. The SW changes of phonons 
consisting in a vibration of intra-trilayer ions or ions in the outer Cu02 planes 
can be elucidated in a similar way. In order to understand the SW changes 
for a composite vibration of the planar oxygens, we have to take into account 
changes of the absolute values of the fields, changes of their relative signs, 
and the displacement pattern. In the frequency region ujq < u p , the local field 
E^ c changes sign under the superconducting transition while the sign of E l °^ 
remains unchanged, i.e., in the SCS the local fields are antiparallel in contrast 
to the NS, where they are parallel. Considering the out-of-phase resonance, 
the eigenvector pattern is in disagreement with that of the local field in the 
NS while they are in accord in the SCS. Consequently, the amplitude of the 
vibration is enhanced under the superconducting transition as well as the SW 
of the phonon [cf. Fig. 8a)]. For the IPR the situation is reversed, i.e., the 
two patterns are in agreement in the NS while they disagree in the SCS. As 
a consequence, the SW of the phonon resonance decreases [cf. Fig. 8 c)]. For 
ujq > u p , the two local fields are parallel in both NS and SCS, and upon 
entering the SCS, their amplitudes increase. Since the eigenvector pattern of 
the in-phase resonance agrees with that of the electric field, its SW increases 
[cf. Fig. 8d)]. 

The shifts of the phonon frequencies are related to the depolarization fields. 
Since the local fields shown in Fig. 9 have been obtained within a model, that 
does not incorporate the phonons, i.e., does not contain the depolarization 
fields, this picture cannot explain the shifts. 



3 Comparison with experimental data and discussion 

Pronounced phonon anomalies have been observed in the far-infrared spectra 
of the trilayer compound Tl-2223 [7]. In this section we discuss the experi- 
mental data, present an assignment of the phonon modes, and compare the 
experimental results with our predictions. Finally we comment on the phonon 
anomalies in Bi 2 Sr 2 Ca2Cu30io (Bi-2223). 

3.1 Tl 2 Ba2Ca 2 Cu 3 Oio 

The conductivity spectra of Tl-2223 shown in Fig. 10 contain a structure sit- 
uated on a high conductivity background. The spectra were obtained [7] by 
Kramers-Kronig transformation of reflectivity data of a ceramic sample con- 
sisting of microcrystals with random orientations. It is known, that sharper 
structures in the data of polycrystalline samples of high-T c compounds cor- 
responds largely to the c-axis polarized phonons and that the metallic back- 
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ground originates predominantly from the a-b conductivity (see [1] and Ref. 
therein). The absence of a-b polarized phonons can be seen, e.g., by comparing 
the reflectivity spectra of a ceramic Y-123 ceramic sample [25] with those of 
a c-axis oriented monocrystal [2]. The phonon structure (see Fig. 10) exhibits 
several changes under the superconducting transition: the SW of the phonon 
at 580 cm -1 decreases, a wide additional absorption band appears in the fre- 
quency region around 500 cm -1 , the SW of the phonon at 370 cm -1 increases, 
the frequency of the phonon at 305 cm -1 decreases by about 20 cm -1 , the SWs 
of the four low frequency phonons increase. 

The lattice dynamical study of Tl-2223 [26] predicts eight infrared-active c- 
axes polarized phonons with eigenvectors shown in Fig. 11. Only seven phonon 
resonances, however, can be found in the normal state spectra of Fig. 10 a). 
We suggest that the phonon, that is not seen in the data, is the one with 
the predicted TO frequency of 443 cm™ 1 because we do not find any phonon 
around 443 cm -1 in the normal state spectra and because this mode should 
have a rather low oscillator strength. Having identified the missing phonon, we 
base the assignment of the phonons on the order of the TO frequencies. The 
differences between the measured and the predicted TO frequencies are within 
20% which can be considered as a good agreement. Note however that also the 
results of the shell model calculations should be taken with some precaution. 
For example, the calculations do not consider free charge carriers. 

As discussed in Sec. 2, the temperature evolution of a phonon resonance de- 
pends on whether its frequency is higher or lower than that of the t-JPR. 
We shall therefore identify the t-JPR in the conductivity spectrum and sub- 
sequently explain the observed changes of the phonons. We suggest that the 
additional absorption band appearing in the SCS around 500 cm -1 corresponds 
to the t-JPR. The band is rather wide (more than 100cm -1 ), similarly as in 
underdoped Y-123, and, in addition, located in the frequency region of no 
phonon related feature in the normal state data. 

The major part of the polarization of the phonon at 580 cm -1 (the 583 cm -1 
mode in Fig. 11) is due to the oxygens in the Tl-0 layers. Vibrations of inter- 
trilayer ions have been described in Sec. 2.2, the corresponding model spectrum 
with the appropriate value of the phonon frequency is shown in Fig. 4 b). 
The calculation reproduces the experimentally observed decrease of the SW. 
Note, however, that the absolute values of the calculated and the measured 
conductivities cannot be compared because of the polycrystalline structure 
of the sample. In addition, in the model spectra of the SCS the phonon line 
appears as a shoulder on the background of the t-JPR; furthermore a week 
broad band emerges at even higher frequencies. Both features seem to be 
present also in the experimental data. The latter can be interpreted as a part 
of the t-JPR, which is separated from the main part by the phonon. This 
feature is even more pronounced in the reflectivity data of the bilayer system 
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Tl2Ba2CaCu2 08 [27,28]. In the low temperature spectrum of this compound 
the t-JPR is centered at about 550 cm" 1 but a part of it appears also around 
700 cm -1 , above the phonon at 600 cm -1 . 

The major part of the polarization of the 370 cm -1 mode of Fig. 10 is due to 
an out-of-phase vibration of the oxygens of the outer and the inner copper- 
oxygen planes (see the 383cm -1 mode in Fig. 11). In section 2.3 we have 
described vibrations of the oxygens in the outer Cu0 2 planes, the vibrations 
of those of the inner plane can be studied along the lines of Sec. 2.2. The 
corresponding model spectra are shown in Fig. 6 a) and Fig. 4 c), respectively. 
Both calculations yield a decrease of the SW and a softening of the phonon, 
in contrast with the experimental observation. However, for composite (out- 
of-phase) vibrations our model yields a SW increase without any frequency 
change [see Fig. 8 a)], in agreement with the experimental data. This agree- 
ment supports the local field picture discussed in the end of Sec. 2.4. In par- 
ticular, it demonstrates dramatic spatial variations of the local electric field, 
even its sign changes on the scale of the inter-plane distance. 

The polarization pattern corresponding to the phonon at 305 cm -1 (292 cm -1 
mode in Fig. 11) is more complicated. The major part of its polarization is due 
to vibrations of the Ca ions and the apical oxygens. In section 2.2, we have 
described vibrations of intra-trilayer and inter-trilayer ions, the corresponding 
model spectra are shown in Fig. 4 c) and Fig. 4 a), respectively. Both exhibit 
a softening of the phonon, the former one shows in addition a SW decrease, 
similar to the one occurring in the experimental data. This suggests that 
the contribution of the Ca ions to the mode is more pronounced. The model 
spectrum of Fig. 4 c) shows a softening of about 40 cm" 1 , much larger than 
the experimental value of 20 cm -1 . It is possible that the presence of other 
phonons in a more complete model would yield a lower value of the frequency 
shift. The SW predicted by the model is smaller than the experimental one, 
both in the NS and SCS. A better agreement could be probably obtained by 
considering a more complex polarization diagram. We have also to keep in 
mind the possibility that the mode at 305 cm -1 corresponds to the in-phase 
mode of the planar oxygens (the 443 cm -1 mode in Fig. 11), i.e., that the order 
of the in-phase and the out-of-phase mode is different from what is suggested 
by the shell model calculations. This seems to be the case in Bi-2223 [6]. 
The corresponding model spectrum shown in Fig. 8 b) exhibits a remarkable 
softening of the phonon resonance and the SW decrease as well. 

Since the ceramic sample consist of randomly oriented microcrystals of the 
highly anisotropic material, the interpretation of the data is not straightfor- 
ward. It is important to find out whether the above discussed spectral changes 
are intrinsic features of the crystal or artefacts induced by the polycrystalline 
structure. The reflectivity of polycrystalline Tl 2 Ba2Cu06 with one Cu02 layer 
per formula unit (see Fig. 12 in Ref. [27]) exhibits only very small differences 
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between above and below T c , but the structures over 250 cm -1 in the data of 
the double layer compound (T^E^CaC^Og) and the triple layer compound 
(Tl-2223) change considerably (see Fig. 13 and 15 of Ref. [27]). We are there- 
fore convinced that the main changes of the phonon structures are intrinsic 
effects rather than artifacts. However, some spectral characteristics must be 
influenced by the polycrystalline structure. We suggest that it is the reason of 
the unexpectedly large values of the SW of phonons [7] and the anomalies of 
the three low frequency modes. 



3.2 A remark on the phonon anomalies in BiaS^CavCuzOio 

Recently A. V. Boris et al. [6] reported the far-infrared c-axis spectra of the 
trilayer compound Bi-2223, that has the same crystal structure as Tl-2223. It 
is interesting to find out whether the phonon anomalies of this compound can 
also be explained by using the model introduced above, especially those of the 
two modes that are specific to the trilayer compounds: the in-phase and the 
out-of phase oxygen bond bending modes. A qualitative explanation of the 
SW changes has already been presented in Ref. [6]. It is further interesting to 
explore the differences between the spectra of the two compounds. 

In the SCS spectra of Bi-2223, a broad additional absorption band appears 
around 500 cm" 1 (see Fig. 1 of Ref. [6]) which has been attributed to the t-JPR. 
Results of the shell model calculations suggest that the spectral structures at 
360 cm -1 and at 400 cur 1 correspond to the in-phase and the out-of-phase 
oxygen bond-bending modes, respectively. The SW of the former decreases in 
the SCS, in contrast to the SW of the latter, that increases below T c . The two 
features are in qualitative agreement with our predictions [see Figs. 8 c) and 
8 a), respectively]. The observed softening of the phonon at 360 cm" 1 (about 
5 cm -1 , the value of the shift is doping dependent), however, is considerably 
smaller than in our model calculations. It may be due to a coupling between the 
phonon and the pronounced mode at 305 cm -1 , mediated by the depolarization 
fields, that is not included in our model. 

The main differences between the spectra of Tl-2223 and Bi-2223 are: (a) The 
spectra of Tl-2223 exhibit only three phonon related structures at frequencies 
above 250 cm" 1 . This is in contrast with Bi-2223, where four main phonon 
bands have been found, (b) The mode at 305 cm" 1 in Bi-2223 has a large SW 
and its contribution to the spectra is almost temperature independent. No 
such mode exists in Tl-2223. (c) The spectral weight of the t-JPR in Tl-2223 
seems to be much larger than in Bi-2223 as described below. 

Measurements on a single-crystal allowed the authors of Ref. [6] to obtain 
the SW of the t-JPR as the difference between the integrated SCS and NS 
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spectra. For a Bi-2223 sample with T c = 102 K that has approximately the 
same frequency of the t-JPR as the Tl-2223 sample studied in [7], the SW of 
the t-JPR is approximately 1100 Q^ 1 cm~ 2 . Concerning Tl-2223, we estimate 
that the SW is about 8000 cm" 2 from the area between the 115 K and the 
25 K spectra as shown in Fig. 10 d). The SW in Tl-2223 is thus several times 
larger than in Bi-2223. It must be of course affected by the ceramic nature of 
the sample, nevertheless, it can be seen even by comparison with the strongest 
phonons in the spectrum that the SW is fairly large. The SW of the t-JPR 
predicted by our model is 18 000 f^cm" 2 . This prediction is about twice as 
much as the observed value for Tl-2223, but the value for Bi-2223 is an order 
of magnitude lower. We emphasize, that the SW of the t-JPR in our model is 
not an arbitrary parameter but is determined by the frequency of the t-JPR 
(see Fig. 2). The reason why the infrared spectra of Tl-2223 and of Bi-2223 
differ so much is not yet understood and remains as an open question for 
future investigations. 



4 Summary and conclusions 

We have extended the model of Ref. [14] that accounts for phonon anoma- 
lies occurring in bilayer compounds to trilayer high-T c superconductors and 
discussed several possible phonon anomalies that can occur in these systems. 
The model explains the experimental infrared data of the trilayer high-T c 
compound Tl-2223 obtained by Zetterer et al. [7]. The broad maximum that 
appears in the conductivity spectra below T c and represents the dominant 
feature of the low temperature spectra corresponds to the transverse Joseph- 
son plasma resonance. The phonon anomalies can be explained in terms of 
changes of the local electric fields induced by the formation of the resonance. 
Of particular interest is the anomaly of the 370 cm -1 mode, i.e., the large in- 
crease of its spectral weight below T c . We recall that this mode consist in an 
out-of-phase vibration of the planar oxygens and is thus specific to the trilayer 
compounds. Its anomalous behavior is due to the fact that the displacement 
pattern does not agree with the pattern of the local field in the normal state, 
but it does agree with the one of the superconducting state, where the sign of 
the field changes from the inner Cu02 plane to the outer one. 

The qualitative and in some cases even quantitative agreement between the 
predictions of our model and the experimental data indicates that the model 
captures well the underlying mechanism. According to the model, the anoma- 
lies are caused by a change of the electronic ground state associated with a 
crossover of the c-axis dynamics from incoherent to coherent, the change which 
would presumably occur even in the absence of any electron-phonon coupling. 
It seems thus that the anomalous infrared- active phonons are passive rather 
than active players in the mechanism of superconductivity. 
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The spectra of Tl-2223 are qualitatively similar to those of Bi-2223 reported 
by Boris et ai: both exhibit the additional absorption band and the anomalies 
of three high-frequency phonons. The magnitudes of the effects, as far as we 
can conclude based on the data of polycrystalline Tl-2223, are considerably 
different. This suggests that the spacing layers of the two materials possess 
fairly different electronic properties. 
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A Appendix 

Here we present the formulas of the approach where boundaries are associated 
with the charged planes of the planar oxygens. Equations (20) and (21) remain 
as they are with 

_ _ / Vti + Vcu \ 24i + d int 
\ Vo ) 44i 



where rjn (rj-mt) is the total surface charge density of the intra-trilayer (inter- 
trilayer) region, i.e., the sum of the densities of the charged planes located 
inside this region. Further 770 (vo) is the surface charge density of the Cu 
ions (O ions) of a Cu0 2 plane. Equation (22) is substituted by 

= -2iuje oX oi(uj)E^:(uj) (A.3) 



with 



0«> = * + # (A.4) 

^0600 ZCZbl + Clint 



The formulas for a, (3, and E l °^p of a bilayer system can be obtained by 
substituting 2du by du- 
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Figure captions: 



Fig. 1. a) Schematic representation of the model of a trilayer superconductor. The 
CuC>2 planes are approximated by homogeneously charged planes. The symbol k 
(— k) denotes the surface charge density of the upper (lower) CuC>2 plane, b) The 
basic structural element of Tl2Ba2Ca2Cu30io [22]. The oxygens in the outer (inner) 
Cu02 planes are labelled as 01 (04), the oxygens in the BaO (T10) planes as 02 
(03). 

Fig. 2. Real part of the c-axis conductivity of the electronic part of the model for 
the following values of the intra-trilayer plasma frequency cjw. (normal state), 
400, 700, 1000 and 1400 cm -1 . It can be seen that Wbi determines not only the 
frequency of the transverse Josephson plasma resonance, but also its spectral weight. 

Fig. 3. Schematic representation of our approach to treat the influence of the phonon 
polarization on the electric field. The plane of ions of a given type is approximated by 
the uniformly charged plane with the same charge per unit area ij. The displacement 
of the plane from its equilibrium position (dashed line) is denoted by r. The electric 
field due to the charged plane is labelled as ±E V , E v = 77/(260600). 

Fig. 4. Anomalies of a phonon consisting in a vibration of ions situated in the 
inter-trilayer region [a) and b)] or in the intra-trilayer region [c) and d)]. The solid 
(dashed) lines correspond to the normal (superconducting) state. The transverse 
Josephson plasma resonance is labelled as t-JPR and the phonon resonance as P03 
and Pea in Figs, a), b) and c), d), respectively. The values of the parameters used 
in the computations are given in the text. 

Fig. 5. The scheme used for the evaluation of the depolarization fields of the oxy- 
gen-bond bending mode. The oxygens of an outer Cu02 plane are represented by 
two homogeneously charged planes. 

Fig. 6. Anomalies of a phonon consisting in vibrations of the oxygens of the outer 
Cu02 planes. The dashed (solid) lines correspond to the normal (superconducting) 
state. The transverse Josephson plasma resonance is labelled as t-JPR and the 
phonon resonance as Poi- The values of the parameters used in the computations 
are given in the text. 

Fig. 7. String model used to describe the composite vibrations of the planar oxygens. 
The displacement of the oxygens 01 (04) from their equilibrium position is labelled 
as r (z), the string constants are denoted by k\ and fo. The local electric fields E^ 
(£ b <f ) act on the oxygens 01 (04). 
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Fig. 8. Anomalies of the in-phase resonance (IPR) and the out-of-phase resonance 
(OPR) of the planar oxygens. The dashed (solid) lines correspond to the normal 
(superconducting) state. The transverse Josephson plasma resonance is labelled as 
t-JPR. The values of the parameters used in the computations are given in the text. 



Fig. 9. Frequency dependencies of the normalized local electric fields in the super- 
conducting state. The symbols E\°%, E^ c and -E^p denote the local fields in the 
inter-trilayer region, inside the trilayer, and on the outer cuprate planes, respec- 
tively. The normalized local fields acting on the ions participating in the modes at 
305, 380 and 580 cm -1 , are indicated by A, B, and C, respectively. 



Fig. 10. Conductivity spectra of Tl2Ba2Ca2CusOio (T c = 112 K) adapted from 
Ref. [7]. The spectra have been obtained by Kramers-Kronig transformation of re- 
flectivity data of a ceramic sample. Part d) contains the spectra of part b) and c). 
The shaded area between the 115 K and the 25 K spectra gives us an rough estimate 
of the spectral weight of the transverse Josephson plasma resonance. 



Fig. 11. Eigenvector diagrams of the c-axis polarized infrared-active phonons of 
Tl2Ba2Ca2Cu30io with the corresponding TO (LO) frequencies obtained by shell 
model calculations in Ref. [26]. 
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Dubroka et al. Fig. 1 
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Dubroka et al. Fig. 2 
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Dubroka et al. Fig. 3 
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Dubroka et al. Fig. 4 
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Dubroka et al. Fig. 5 
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Dubroka et al. Fig. 6 
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Dubroka et al. Fig. 7 
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Dubroka et al. Fig. 8 
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Dubroka et al. Fig. 9 
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Dubroka et al. Fig. 10 
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Dubroka et al. Fig. 11 






j45^ Jtt Jtt 





J&J? 

©4 © 




o 




k& 




o 

«) (» 



72.8 (103) 104 (131) 150 (151) 



o 






175 (197) 



© 




292 (372) 383 (421) 443 (454) 583 (616) 



32 



